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Abstract 

It is proved that the two-photon Rabi Hamiltonian is quasi exactly 
solvable on the beisis of the two diflterent quadratic Lie algebras. 

1 Introduction 

By definition, a quasi-exactly solvable (QES) Hamiltonian H has a finite- 
dimensional matrix representation within a finite-dimensional subspace of the 
whole infinite-dimensional Hilbert space. As a result, a partial algebraization 
of the spectrum occurs. Correspondingly, a QES Lie algebra is the Lie algebra 
of differential operators g {Ji, J2, Jn) that admits finite-dimensional repre- 
sentation in the subspace 3?Ar of smooth functions { if ip £ SR^v and Jn G 
then J„ (ip) g 5ft AT ) [T]. Then, the QES Hamiltonian belongs to an universal 
enveloping QES algebra. 

In investigation of QES systems, it was implied up to now that such a system 
has only one underlying QES algebra. In the present paper, we would like to 
pay attention that a given QES Hamiltonian may have more than one QES 
algebra. In doing so, different finite-dimensional representations of such QES 
algebras gi, 02, 03.... cut out different parts of the spectrum of Hamiltonian 
H and/or correspond to different pieces in the space of its parameters. This 
observation not only deepens our understanding of the mathematical structure 
of QES systems but is also important for physical applications since it allows to 
explore the properties of the spectrum in more detail. 

In the present paper, we consider the two-photon Rabi Hamiltonian (TPRH). 
We find explicitly two different QES algebras that ensure its quasi-exact solv- 
ability for certain values of parameters. It turns out that both the algebras 
are quadratic in contrast to the case of linear sl(2,C) algebra typical of one- 
dimensional systems [2] - [2] . ( Quadratic algebras in the QES context were also 
discussed in Ref. [10] ). 



The fact that TPRH is QES was shown in Ref. [H]- The works [H], [12] 
and the present one deal with different relationships between the parameters of 
the system and describe, in general, different parts of the spectrum for which 
different QES Lie algebras are relevant. Such a rich QES structure of TPRH is 
a separate interesting property of this system. 

2 Two-photon Rabi Hamiltonian 

Rabi Hamiltonian [13 describes a two-level system (atom) coupled to a single 
mode of radiation via dipole interaction. It reads 

i^R = YT,+Li.6+fe + g(6 + fe+) (1) 

where b and 6"'' are the Bose annihilation and creation operators respectively 
([6,6^] = 1) and a^, a± = ax i i ■ cry are the Pauli matrices. It is often 
considered in the rotating-wave approximation in which the terms proportional 
to 6(T_ and 6^"(t+ are omitted in which case it is sometimes called the Jaynes- 
Cummings Hamiltonian [14j . [15j (in some works, the latter term applies to the 
full Hamiltonian (fT|) \W\). 

The simplest generalization of such a Hamiltonian taking into account non- 
linear optical processes is TPRH: 

H = +[b-b+b + g (b^ + (6+)') • (a+ + a_) (2) 

The Hamiltonian ^ describes the same physical system "two-level atom 
plus radiation" but emission and absorption of radiation occurs due to a two- 
photon process instead of a single-photon ones in the ordinary Rabi Hamiltonian 
([13|V In the present work, we consider the full form of Hamiltonian ^ without 
using the rotating wave approximation. 

In what follows we will use the representation of the Pauli matrices 

= ( -1 ) ' '"^ = ( ) ' '"^ = ( 1 ) ■ 

Then, the spectral problem Hijj = Etjj for Hamiltonian ^ can be represented 
in the form 

where L± = b+b ± g/2 ■ ({b+ f + b^^ , ujq = Cuq/ (2 • w), g = 4 • g/w, E = E/Co. 

Eliminating the function ip2 from the system of equations ([3]) we arrive at 
the fourth-order equation for the component Lpi of the wave function: 

Liipi = 

Li = {L^L+-E{L+ + L^) + E^ ^ul) (4) 
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We find the conditions under which the operator Li is QES and, hence, 
so is Hamiltonian ([2]). We would hke to draw attention to the following sub- 
tlety. We are dealing with the two-component wave function. Usually, while 
considering matrix QES Hamiltonians [T7] - [5D] it was assumed that all com- 
ponents of the wave function belong to the invariant subspace. Now, this is, 
generally speaking, not the case. We require that, say, the component <pi obey 
QES equation. However, the second component ip2 = 1 / i^o ■ {E — L+) ipi does 
not, in general, belong to the invariant subspace. (We may interchange the role 
of ipi and (p2 and require that (p2 obey the QES equation.) 

Now we will show that the operator Li ^ can be represented in terms of 
generators of two QES algebras and, therefore, the operator Li itself possesses 
the same invariant subspaces as the aforementioned generators. The explicit 
description of corresponding generators and their invariant subspaces is given 
in Appendix A. The fact that the generators under discussion are quasi-exactly 
solvable is by no means obvious and was established within the approach called 
in [12] QES-extension. It was shown there that the differential operators J2 
([^ possess invariant subspaces 3?^ , ^% which are not polynomials 
and the generators cannot be obtained as a polynomial deformation of the linear 
algebra sl{2,R) [TU]. It turned out that they can be expressed in terms of 
hypergeometric functions [12] . I use here these results from the previous work 
[12] with the notations borrowed from that paper. It is also interesting that the 
corresponding algebras are quadratic in terms of aforementioned operators as it 
is seen from the commutation relations ((28)) . (|3T]) . 

Let us use the coordinate representation 

b=^,b+ = ^^ (5) 

in which 2 ■ L± = {±g — 1) ■ + (1 ± g) • — 1. We want to express Li ^ in 
terms of operators J2 , J 2 or ^ J3 ■ As these operators have the second order 
in derivatives while Li is of the fourth order, we are seeking for the suitable 
quadratic combinations. Using also the gauge freedom, we rely on the relation 

e/3(-).L,.e-/'(-)=P(j-,J+)|,^,(,) (6) 

(similarly for J3" , J^) where (3{x) is a gauge function, P {z, y) ^ c^z ■ +Cyy ■ 
+Cyz • yz +Czy • zy +Cz ■ z +Cy ■ y +cq — is the second order polynomial, 77 (x) 
is the function that defines the change of variable. Making a simple choice 
P (x) — const ■ x^, r] (x) = const ■ x^ we obtain the system of equations for the 
coefficients of the polynomial P {z,y). Below we give several solutions which we 
managed to obtain. 

I) The subspace 3?^: 



Ag'-{j-f+4g-S2^2g-'-J-+iJ+-u;-'o-ao (7) 



where ap = (3 + 47V) • (47V {g^ - 1} - S + bg^) /4, 

a = —3/4, s = —1/2 ( the parameters of subspace ) (8) 
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the parameter of change of variable ) 



(9) 
(10) 



E = -^ + 2{N +1) ■ \/Y~g^ ( the energy of the Hamiltonian © ). (11) 



II) The subspace 3?^ 



N 



43' • ( J3 



4 Jg*" — cJq — ao 



where ao = (1 + 2iV) • (2iV (5^ _ 1) _ 1 + 3^2) 
3 

a = - — ^,s=:l/2( the parameters of subspaces 3?' 



1 1-5 



2V 1+5 
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( the parameter of gauge transformation ) 
: ( the parameter of change of variable ) 



(12) 

(13) 

(14) 
(15) 



E=-- + {N+l)- Vl - 5^ ( the energy of the Hamiltonian ([3]) ). (16) 

The expressions ([7]), ([T2|) demonstrate explicitly that Hamiltonian ^ is 
indeed QES for corresponding values of parameters. Using the asymptotic for- 
mulas for the hypergeometric function [21], one can check that for the values of 
parameters ([9]), pO)) . (fT4|) . ^5]) the functions belong to the space L2 (—00, +00). 
For example, for the subspace 3?^ for a; ^ cxd we have (fi ^ exp (— (c — ^) a;^) 
[H] where c-C>0, 0<g<l(II|. 

Now, the matrix representation of the operator Li ^ follows from those of 

(PS)) and jf- (1^^ . Below we list explicit formulas for small values of iV . For 
convenience, we also include in the subspaces , 3?^, the gauge factor e~'^^ 
fromeqs. ([7|). (fH]) . 

I) The subspace (TV = 0): 



Q — span < e 



iFi 



-1/2 ' ''■^ 



1/2' "^-^ 



(17) 



where c and ^ are the parameters pl7|) . Finite-dimensional representations 
of the operators J>^ take the form: 



1/0 

1 -1 



1 / -3 6 
1 



(18) 
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whence, taking into account eq. ([7]) we obtain the matrix representation for the 
operator Li on the subspace : 



Li 



3.9 + 9/4 - cj2 



-35(25 + 1) 



2{g + l) 1/4 - 35 



(19) 



The spectral problem ^ possesses a non-zero solution if the determinant of the 
matrix (fTO|) is equal to zero. The only real solution g (ojq) of the corresponding 
equation that obeys the condition < g < 1 has the form 



' (^o) 



V(4c^g-9)(1 
8^0 



4^2) 1 3 



(20) 



The components of the wave function corresponding to the solution ([20)1 for the 
particular values of parameters are listed in Appendix B. 
II) The subspace ^ {N ^ 2): 



iFi 



span • 



1/2 '"^-^ 



1/2'^-^ 



1/2'^-^ 



(21) 



where c and are the parameters (|T4l) . (fT5)l . Taking into account ((30)) . p2| we 
find the finite-dimensional representations of operators J^: 



(22) 



whence, with l|12p taken into account, we also obtain the matrix representation 
for the operator Li on the subspace 3?2 : 






2(9 + 1) 

33/4 - Loi - 8g^ -3 [g - 
10(l-<?) -15/4- 



(23) 



The spectral problem ([4]) has a non-zero solution if the determinant of the matrix 
(|23|) is equal to zero. We find two solutions g± (wq) satisfying the condition 
< 5± ('^o) < 1: 



g± (wo) = z 



1 /(5t2wo)(2wo±3)(1±2wo) 1 



1 



-<wo<2±-. (24) 



The components of the wave function corresponding to the solution 9+ (o^q) are 
listed in Appendix B for particular values of parameters. 

Equations (O, obtained above give evidence that Hamiltonian ^ is 
quasi-exactly solvable on the basis on two quadratic Lie algebras (see Appendix 
A). Each of them cuts out some set of eigenvalues from the spectrum of the total 
Hamiltonian. Thus, different parts of the spectrum are described by different 
Lie algebras. 
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3 Summary and outlook 



It is shown that the two-photon Rabi Haniihonian ([2]) is quasi-exactly solvable 
on the basis of two different quadratic Lie algebras. The results of the present 
paper demonstrate that different parts of the spectrum correspond to different 
Lie algebras. It is especially interesting that the functions which realize invariant 
subspaces are expressed in terns of special functions (cf. also [2]) whereas 
application of the linear Lie algebra leads to subspaces spanned on polynomials 

m- 

The results obtained provoke also further questions which we outline here 
briefly. First of all, it concerns the comparison of the present article with previ- 
ous results. Both in Ref. [TT] and in our paper QES solutions are obtained for 
the energies 

E = -^ + M-^l-g^ (25) 

where M is a half- integer in JTj and an integer in our paper. In [11], these 
energies were claimed to correspond to level-crossing of the TPRH but not to 
all of them. It was conjectured there that integer M might correspond to the 
remaining crossings. Now, in view of the present results, it would be especially 
interesting to elucidate the relationship between level-crossings and different 
types of quasi-exact solvability. 

In the previous paper |12j we used another representation - the Fock-Bargmann 
one instead of the coordinate one in the present work. Correspondingly, the op- 
erators Li and its analogue L from [1^ do not coincide, so that the conditions 
of their quasi-exact solvability impose different restrictions on the parameters 
of Hamiltonian. It turned out that in the first case it reduces to g* = ^[12] 
independently of N whereas in our case g is the function of ujq, the concrete 
form of this function being different for different N (see eq. (PH)) for iV = and 
(gl for = 2). 

It is of separate interest the issue of classification of states. In particular, 
it concerns the question as to which multiplet the ground state belongs. If 
the corresponding relevant QES algebra is known, the problem of minimiza- 
tion of the functional (^^,H'i!^ reduces from the whole Hilbert space to its 
finite-dimensional part that simplifies greatly the problem and can be useful for 
applications. 
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5 Appendix A. 

Here we give explicit formulas for two QES quadratic Lie algebras used in the 
text [H]. 
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I) Finite dimensional function subspace TZj^ — span{fQ , /j^, /q , 
/-} (TV = 0,1,2,....), dim{n%) = 2{N+l), formed by functions f+ = ■ 

1^1 (n = 0,1,.., TV -1, TV) [H], is invariant for 



the operators J^: 

'^2-=i^ + (l + «-0i 
The operators act on the functions f^{x) G 7?.^ as follows: 



(26) 



' \fnj \{S - n) (1 - AO • fn + S„ • + S (2S„ - 1) • /rt/ 



^2" 



/n+\ ^ A" - ") /rt + « + s) • /+_! + a (1 + 2n) /s ■ /„ 
/«/ V (a + n+l)/,7+n(n-s) •/-_! + 2ns 



where A„ = n — 2A^, i?„ = n — N. The commutation relations of the operators 
are: 

[-^2^; -'^2 ]= ■5*2 

[ J2+ , ^2] = 4 • J2+ - 2 • ^2 - 4 J+ - c+ - c+ (28) 
[J2-, ^2] = -2 • (J2-)' - J2+ + c+J.^ ^(a~N){s + 1) 

where c+ = 2 (1 + a)+s, 4 = (2TV - s) (s - 2 - 2TV), 4 = {N - a) (2 + 2iV - s) ■ 
(s + 1). 

II) Finite dimensional function subspace TZf^ — span{fQ, /i, /at} (A^ = 
0, 1, 2, ....), dim {n%) = TV + 1, formed by functions /„ = iFi ; i] (n = 
0, 1, .., TV — 1, iV) j21j, is invariant for the operators J^: 

\n , (29) 

The operators jf' act on the functions /„(a;) G 7?.^ as follows: 

^3^ {fn) = {n + a)- /„ (30) 



•^3^ {fn) = {s-n + Cn) ■ fn + «„ • B„ • /„+i + n (a„ - s) • /„- 



1 



where a„ = n+a, C„ = a„ • (A^ — 2n), _B„ = n— TV . The commutation relations 
of the operators jf are: 

[j3+,53] = 4- J3+. J3--2.53-C+- J2+ + 4-'^2"+C^ (31) 

[J3-, 53] = -2 • {J^f -j++c+.J--s-a 
where = .s + 2a + TV , = (TV + 2 - s) • (TV - s) , c| = sa- (TV + 2 - s). 
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The operators (a, s, N) can be transformed into the operators 
(a', s', N') (|29p by means of substitution of the parameters s s' — 1 , a ^ 
a' + {N' - 1) /2, TV ^ (iV' - 1) /2 

{a' + [N' - 1) /2, s' - 1, {N' - 1) /2) ^ {a', s', N') . (32) 

This substitution has meaning for odd N' only. It is also worth noting that the 
subspace is the subset of the subspace ( C di% ) . This claim follows 
from the properties of the hypergeometric functions i2T] . 



6 Appendix B. 

I) The subspace {N = 0): 

According to eq. ([4]), we need the eigenvector corresponding to the zero 
eigenvalue of the transposed matrix (jl9p where we should take into account 
(|^D|) . Then, we obtain for the wave function: 



* = 



V2 



Pll ■ hi (x) +P12 ■ h2 (x) 
P21 ■ hi (x) +P22 ■ h2 (x) 



(33) 



where hi (x) 



4^, Pll 



iFi 
= 250 



'^^'^■£x^ 

-1/2'?-^ 



, h2 (x) = e-=- i^^i 
68V15, pi2 = -235 - 6OVT5, P21 = 70 + 28^/15, 



2{VT5-l) 

P22 = 35— (120 + 20^/15) x^. The second component <p2 was found from relation 



1 , 

ip2 = — [E - L+) ipi 



(34) 



The solution (|33p (not normalized) corresponds to the energy value E = 5/4. 
II) The subspace 3?^ (N = 2) 

Proceeding along the same lines as in I) and taking into account 
obtain the wave function corresponding to the solution (ujq) at = 1 



we 



* = 



where hi (x) 



f2 

iFi 



Pll ■ hi (x) +P12 ■ h2 (x) 
P21 ■ hi (x) +P22 ■ h2 (x) 



(35) 



3/2 ' ''■^ 



i-t'i 



1/2 '"5"^ 



\/8-3V5 ^ 
2-\/8+3vf' 



, h2 (x) 

3x2 ^7^^o + 3184^/5) x^ - ^19 (371 + 170^5)), 
19 (6368 + 2844 V5) x2 - (l235 + 418%/5)), p2i = f (79 + 32^5) 



-^\/ -Pll 



Pl2 = -- 

P22 (x) = (76x2 — -y/l9 (70 — 31\/5)) . rp]_^g second component ip2 was ob- 
tained from relation ((34| . Solution p5|) (not normalized) corresponds to the 
energy value E = 3\/T9/8 - 1/2. 
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